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Abstract. Large volume cosmological simulations succeed in reproducing the large-scale struc- 

QQ " ture of the Universe. However, they lack resolution and may not take into account all relevant 

f^ I physical processes to test if the detail properties of galaxies can be explained by the CDM 

f~*) ' paradigm. On the other hand, galaxy-scale simulations could resolve this in a robust way but 

fS) do not usually include a realistic cosmological context. 

,_^ To study galaxy evolution in cosmological context, we use a new method that consists in 

^ coupling cosmological simulations and galactic scale simulations. For this, we record merger and 

' , gas accretion histories from cosmological simulations and re-simulate at very high resolution the 

C^ • evolution of baryons and dark matter within the virial radius of a target galaxy. This allows us 

^~H I for example to better take into account gas evolution and associated star formation, to finely 

. study the internal evolution of galaxies and their disks in a realistic cosmological context. 

We aim at obtaining a statistical view on galaxy evolution from z ~ 2 to 0, and we present 

Qj. here the first results of the study: we mainly stress the importance of taking into account gas 

I ■ accretion along filaments to understand galaxy evolution. 
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1. Introduction 



The morphology of galaxies in the Local Universe is well constrained by observa- 
tions, but is still largely unexplained. Indeed, large volume cosmological simulations 
^O fail to reproduce realistic galaxies. For instance, the disks formed are often too con- 

T-H centrated : it is the "angular momentum problem" , well known since the early work 

[^ , of Navarro & Benz (1991) It is still unclear whether this is an intrinsic problem of the 
f~^ ACDM paradigm or if something (i.e. resolution, physical processes...) is missing in these 

OO , simulations. 
^^ ■ Another puzzle is the question of disk survival till z=0 ( |Koda et al., 2007| . For in- 



stance, Kautsch et al. (2006) study a large sample of edge-on spiral galaxies in the SDSS 
and find that a significant fraction of them (i.e. roughly one third) are bulgeless or 
^ , "superthin". This is still unexplained by cosmological models. Indeed, ACDM predicts 
C^ • that galaxy interactions are frequent (see e.g. the recent work by [Stewart et al., 2007D . 
More exactly, major mergers, that are well known to destroy disks to form ellipticals 



(Barnes & Hernquist, 1991) are rather rare, but minor mergers are much more com- 
mon. These minor mergers can thicken disks, and if frequent enough could even form 
elliptical galaxies (Bournaud et al., 2007) . The problem is then to find whether ACDM 
predicts too many mergers, or if the satellites have properties and orbital parameters such 
that they have little influence on the galactic disks. Also, gas accretion along filaments 



could fuel a thin disk and counteract the effect of mergers (Dekel & Birnboim, 2005 



Keres et al., 2005 Ocvirk et al., 2008) 
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To study the properties of galaxies at low and high redshift, it thus seems necessary 
to take the full cosmological context into account. Large scale cosmological simulations 
could of course achieve this goal and give a statistical view on galaxies at each redshift, 
but for now they mainly lack resolution at the galactic scale. On the contrary, small 



volume cosmological simulations like the one performed by Naab et al. (2007) can resolve 
galactic scales in detail but are so time-consuming that obtaining a statistical sample is 
for now a challenge. 

A first method to solve these problem is to use semi-analytical models, i.e. extracting 
merger trees from cosmological simulations and using different recipes to infer physical 



properties of galaxies ( Somerville et al., 2001 Hatton et al., 2003 Khochfar fc Burkert, 2005 ) 



The drawback is that approximations are necessary. 

Another possibility has been explored by Kazantzidis et al. (2007)[ Read et al. (2007) 



and Villalobos & Helmi (2008) : they extract merger histories from cosmological simu- 
lations and re-simulate these histories at higher resolution. Nevertheless, they perform 
collisionless simulations with no gas component, neither in the main galaxy, nor in satel- 
lites, nor in filaments. 

We here present a new approach where we re-simulate at high resolution a history given 
by a cosmological simulation, using self consistent realistic galaxies (the main galaxy and 
the satellites have a gas disk, a stellar disk and a dark matter halo), and we also take into 
account gas accretion from cosmic filaments. Our goal is to obtain a statistical sample 
of merger and accretion histories in a ACDM context to simulate the resulting galaxies 
and to compare our results to observations at various redshifts. 

After a description of the technique used, we will present our first results and emphasize 
the importance of gas accretion along filaments to understand galaxy evolution. 



2. Method 

2.1. Analysis of the cosmological simulation 

Merger histories and accretion data are extracted from a dark matter only cosmological 
simulation performed with the AMR code RAMSES ( [Teyssier, 2002p . This simulation 
has an effective resolution of 512'^ and a comoving box length of 20 h~^ Mpc. The 
mass resolution is 6.9x10® M©, so that a Milky Way type halo is made of a few 10^ 
particles. The cosmology is set to ACDM with J7,„=0.3, fl\—0.7, Ho=70 km.s~^.Mpc~^ 
and (78=0.9. 



In this simulation, halos are detected with the HOP algorithm (Eisenstein & Hut, 1998), 
with (5pcak=240, 4addio=200 and (5outoi=80 (the minimal number of particles per halo is 
fixed to 10). In the following, we also take into account particles that do not belong to a 
halo, and we consider them as diffuse accretion. 

The halo of which we want to build the merger and accretion history is then chosen 
in the final snapshot of the simulation (at z = 0) and is traced back to higher redshift 
(typically z ~ 2) : we will call it the main halo. From z ~ 2 to z = 0, each halo or 
particle (in the case of diffuse accretion) entering a sphere around the main halo (the 
radius of this sphere is the virial radius of the main halo at z=0) is recorded, with its 
mass, position, velocity and spin (spin is of course omitted for diffuse accretion). 

2.2. High resolution re-simulation 
2.2.1. The PM code 

The history that has been extracted from the cosmological simulation is re-simulated 



with a particle-mesh code (Bournaud & Combes, 2002). Gas dynamics is modeled with a 
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Figure 1. Left : Initial distribution of stars (top panel) and gas (bottom panel) for the main 
galaxy, seen face-on and edge-on (each panel is 40 kpc x 40 kpc in size) . Right : large scale view 
of the gas distribution in a simulation box (the panel is 440 kpc x 440 kpc in size) . 

sticky-particle scheme with f3r=0-8 and Pt—0-7, and star formation is computed according 
to a Kennicutt law with an exponent 1.5. 

The maximum spatial resolution is 130 pc. For the two simulations shown hereafter, 
the mass resolution varies from 1.2x10^ Mq to 2.1x10'' M© for gas particles, from 6x10"^ 
M0 to 1.4x10^ Mq for stellar particles and from 1.2x10^ Mg to 4.4x10^ Mq for dark 
matter particles. This allows to have a total number of particles of the order of 15x10^ 
at the end of both simulations. 



2.2.2. Model galaxies 

Each halo of the cosmological simulation (i.e. the main halo as well as all the interacting 
satellites) is replaced with a realistic galaxy, having a disk, a bulge and of course a dark 
matter halo. The total mass of the galaxy is divided in 20% of baryons and 80% of dark 
matter (the mass of dark matter being given by the cosmological simulation) . The dark 
matter halo follows a Burkert profile extended to its virial radius, with a core radius 



chosen to follow the scaling relations given in Salucci & Burkert (2000) The disk radius 



of each galaxy is proportional to the square root of its mass so that the surface density 
is constant from one galaxy to another. 

The gas fraction in the disk is 30% for galaxies that have a halo mass lower than 
10^^ Mg. For galaxies that have a greater halo mass, the gas fraction is set to 30% at 
high redshift (z>0.8) and 15% at low redshift. 

Figure [1] (left side) shows for example the initial distribution of gas and stars in the 
main galaxy. 



2.2.3. Diffuse accretion 

Each dark matter particle that is considered as diffuse accretion in the cosmological 
simulation is replaced with a small blob of particles, containing in mass 20% of gas and 
80% of dark matter. 

The right side of figure [1] shows an example of simulation where the main galaxy (edge- 
on) is surrounded by accreted gas (clearly in a filament) and a few satellite galaxies. 
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Figure 2. Evolution of tlie total mass of dark matter in the simulation box as a function of 
time for the two simulations studied here : in the left case, the mass growth is dominated by 
accretion, and in the right one by mergers 

2.3. Two examples 

We present here the first results concerning two simulations, that have been chosen to 
have a mass at z=0 of the order of magnitude of the mass of the Milky Way. They have 
very different histories. 

In the the first one, the mass growth of the galaxy is dominated by diffuse accretion 
(at a mean rate of ~ 5 Mq/jv). Only some very minor mergers take place, the most 
important of these mergers having a mass ratio of 12:1 (see on the left panel of figure [2] 
the mass evolution as a function of time). We will call this simulation "the calm case". 

The other simulation also contains diffuse accretion, but is mainly dominated by merg- 
ers. There is a first period of repeated minor and major mergers (mass ratios 8:1, 10:1, 
3:1 and 4:1) at the very beginning of the simulation, then a calm phase and finally a 
major merger (mass ratio 1.5:1) at low redshift (see right panel of figure[2]). We will call 
it "the violent case". 



3. Results 

3.1. The calm, case 

The evolution of the distribution of gas and stars is shown in figure [3l Gas is smoothly 
accreted around the galaxy and falls onto the disk. Minor mergers are not strong and 
frequent enough to destroy the stellar disk. They only slightly heat it, and a thin stellar 
disk is rebuilt thanks to gas from diffuse accretion along the filaments. The thin disk is 
mainly formed from stars younger than 4 Gyr, and has a well-defined structure with two 
spiral arms. 

3.2. The violent case 

In this case, the evolution of the morphology of the galaxy is totally different (see figure 
13]). The disk is destroyed early by the first series of mergers. In fact, after the first of 
these mergers (which has a mass ratio of 8:1) the disk is already very perturbed, and the 
following mergers contribute to the transformation of the galaxy into an elliptical. 

Nevertheless, thanks to gas accretion that takes place along a filament, a gas disk is 
gradually re-built into the elliptical galaxy (this would not happen if only mergers were 
taken into account in the simulation). New stars form in this disk, forming a young stellar 
disk inside the old spheroid (see figure [5]) , this disk being in a perpendicular plane with 
respect to the initial disk. Finally, the last major merger (with a mass ratio of 1.5:1) 
destroys this disk and the galaxy becomes elliptical again. 
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Figure 3. Evolution of tlie distribution of gas (top panels) and stars (bottom panels) for the 
calm case. Snapshots are taken every Gyr and each panel is 40 kpc x 40 kpc in size. 

4. Conclusion 

In order to study galaxy evolution in cosmological context, we have successfully de- 
veloped a technique that allows us to perform high resolution simulations taking into 
account realistic merger and gas accretion histories. 

The first two simulations shown here do not allow us to draw any general conclusion 
on galaxy evolution in a ACDM context. Nevertheless, we can already confirm that even 
low mass satellites can thicken disks and that ellipticals from both through repeated 
minor mergers and major mergers. We also emphasize that gas accretion from filaments 
can allow to rebuild a thin disk in a galaxy, which proves the absolute necessity to take 
this accretion into account to understand galaxy evolution. 
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Figure 4. Evolution of tlie distribution of gas (top panels) and stars (bottom panels) for the 
violent case. Snapshots are taken every Gyr and each panel is 40 kpc x 40 kpc in size. 




Figure 5. Projected stellar mass density at z = 0.2 for the violent case. 
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